3.380 \(\int \frac{\sqrt{1-x^2}}{x^3 (a+b x^2+c x^4)} \, dx\)

Optimal. Leaf size=290 \[ -\frac{\sqrt{c} \left (\frac{a (b-2 c)+b^2}{\sqrt{b^2-4 a c}}+a+b\right ) \tanh ^{-1}\left (\frac{\sqrt{2} \sqrt{c} \sqrt{1-x^2}}{\sqrt{-\sqrt{b^2-4 a c}+b+2 c}}\right )}{\sqrt{2} a^2 \sqrt{-\sqrt{b^2-4 a c}+b+2 c}}-\frac{\sqrt{c} \left (-\frac{a (b-2 c)+b^2}{\sqrt{b^2-4 a c}}+a+b\right ) \tanh ^{-1}\left (\frac{\sqrt{2} \sqrt{c} \sqrt{1-x^2}}{\sqrt{\sqrt{b^2-4 a c}+b+2 c}}\right )}{\sqrt{2} a^2 \sqrt{\sqrt{b^2-4 a c}+b+2 c}}+\frac{(a+2 b) \tanh ^{-1}\left (\sqrt{1-x^2}\right )}{2 a^2}-\frac{1}{4 a \left (1-\sqrt{1-x^2}\right )}+\frac{1}{4 a \left (\sqrt{1-x^2}+1\right )} \]

[Out]

-1/(4*a*(1 - Sqrt[1 - x^2])) + 1/(4*a*(1 + Sqrt[1 - x^2])) + ((a + 2*b)*ArcTanh[Sqrt[1 - x^2]])/(2*a^2) - (Sqr
t[c]*(a + b + (b^2 + a*(b - 2*c))/Sqrt[b^2 - 4*a*c])*ArcTanh[(Sqrt[2]*Sqrt[c]*Sqrt[1 - x^2])/Sqrt[b + 2*c - Sq
rt[b^2 - 4*a*c]]])/(Sqrt[2]*a^2*Sqrt[b + 2*c - Sqrt[b^2 - 4*a*c]]) - (Sqrt[c]*(a + b - (b^2 + a*(b - 2*c))/Sqr
t[b^2 - 4*a*c])*ArcTanh[(Sqrt[2]*Sqrt[c]*Sqrt[1 - x^2])/Sqrt[b + 2*c + Sqrt[b^2 - 4*a*c]]])/(Sqrt[2]*a^2*Sqrt[
b + 2*c + Sqrt[b^2 - 4*a*c]])

________________________________________________________________________________________

Rubi [A]  time = 2.36223, antiderivative size = 290, normalized size of antiderivative = 1., number of steps used = 8, number of rules used = 6, integrand size = 29, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.207, Rules used = {1251, 897, 1287, 207, 1166, 208} \[ -\frac{\sqrt{c} \left (\frac{a (b-2 c)+b^2}{\sqrt{b^2-4 a c}}+a+b\right ) \tanh ^{-1}\left (\frac{\sqrt{2} \sqrt{c} \sqrt{1-x^2}}{\sqrt{-\sqrt{b^2-4 a c}+b+2 c}}\right )}{\sqrt{2} a^2 \sqrt{-\sqrt{b^2-4 a c}+b+2 c}}-\frac{\sqrt{c} \left (-\frac{a (b-2 c)+b^2}{\sqrt{b^2-4 a c}}+a+b\right ) \tanh ^{-1}\left (\frac{\sqrt{2} \sqrt{c} \sqrt{1-x^2}}{\sqrt{\sqrt{b^2-4 a c}+b+2 c}}\right )}{\sqrt{2} a^2 \sqrt{\sqrt{b^2-4 a c}+b+2 c}}+\frac{(a+2 b) \tanh ^{-1}\left (\sqrt{1-x^2}\right )}{2 a^2}-\frac{1}{4 a \left (1-\sqrt{1-x^2}\right )}+\frac{1}{4 a \left (\sqrt{1-x^2}+1\right )} \]

Antiderivative was successfully verified.

[In]

Int[Sqrt[1 - x^2]/(x^3*(a + b*x^2 + c*x^4)),x]

[Out]

-1/(4*a*(1 - Sqrt[1 - x^2])) + 1/(4*a*(1 + Sqrt[1 - x^2])) + ((a + 2*b)*ArcTanh[Sqrt[1 - x^2]])/(2*a^2) - (Sqr
t[c]*(a + b + (b^2 + a*(b - 2*c))/Sqrt[b^2 - 4*a*c])*ArcTanh[(Sqrt[2]*Sqrt[c]*Sqrt[1 - x^2])/Sqrt[b + 2*c - Sq
rt[b^2 - 4*a*c]]])/(Sqrt[2]*a^2*Sqrt[b + 2*c - Sqrt[b^2 - 4*a*c]]) - (Sqrt[c]*(a + b - (b^2 + a*(b - 2*c))/Sqr
t[b^2 - 4*a*c])*ArcTanh[(Sqrt[2]*Sqrt[c]*Sqrt[1 - x^2])/Sqrt[b + 2*c + Sqrt[b^2 - 4*a*c]]])/(Sqrt[2]*a^2*Sqrt[
b + 2*c + Sqrt[b^2 - 4*a*c]])

Rule 1251

Int[(x_)^(m_.)*((d_) + (e_.)*(x_)^2)^(q_.)*((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4)^(p_.), x_Symbol] :> Dist[1/2,
Subst[Int[x^((m - 1)/2)*(d + e*x)^q*(a + b*x + c*x^2)^p, x], x, x^2], x] /; FreeQ[{a, b, c, d, e, p, q}, x] &&
 IntegerQ[(m - 1)/2]

Rule 897

Int[((d_.) + (e_.)*(x_))^(m_)*((f_.) + (g_.)*(x_))^(n_)*((a_.) + (b_.)*(x_) + (c_.)*(x_)^2)^(p_.), x_Symbol] :
> With[{q = Denominator[m]}, Dist[q/e, Subst[Int[x^(q*(m + 1) - 1)*((e*f - d*g)/e + (g*x^q)/e)^n*((c*d^2 - b*d
*e + a*e^2)/e^2 - ((2*c*d - b*e)*x^q)/e^2 + (c*x^(2*q))/e^2)^p, x], x, (d + e*x)^(1/q)], x]] /; FreeQ[{a, b, c
, d, e, f, g}, x] && NeQ[e*f - d*g, 0] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - b*d*e + a*e^2, 0] && IntegersQ[n,
 p] && FractionQ[m]

Rule 1287

Int[(((f_.)*(x_))^(m_.)*((d_) + (e_.)*(x_)^2)^(q_.))/((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4), x_Symbol] :> Int[Ex
pandIntegrand[((f*x)^m*(d + e*x^2)^q)/(a + b*x^2 + c*x^4), x], x] /; FreeQ[{a, b, c, d, e, f, m}, x] && NeQ[b^
2 - 4*a*c, 0] && IntegerQ[q] && IntegerQ[m]

Rule 207

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> -Simp[ArcTanh[(Rt[b, 2]*x)/Rt[-a, 2]]/(Rt[-a, 2]*Rt[b, 2]), x] /;
 FreeQ[{a, b}, x] && NegQ[a/b] && (LtQ[a, 0] || GtQ[b, 0])

Rule 1166

Int[((d_) + (e_.)*(x_)^2)/((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[b^2 - 4*a*c, 2]}, Di
st[e/2 + (2*c*d - b*e)/(2*q), Int[1/(b/2 - q/2 + c*x^2), x], x] + Dist[e/2 - (2*c*d - b*e)/(2*q), Int[1/(b/2 +
 q/2 + c*x^2), x], x]] /; FreeQ[{a, b, c, d, e}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - a*e^2, 0] && PosQ[b^
2 - 4*a*c]

Rule 208

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[-(a/b), 2]*ArcTanh[x/Rt[-(a/b), 2]])/a, x] /; FreeQ[{a,
b}, x] && NegQ[a/b]

Rubi steps

\begin{align*} \int \frac{\sqrt{1-x^2}}{x^3 \left (a+b x^2+c x^4\right )} \, dx &=\frac{1}{2} \operatorname{Subst}\left (\int \frac{\sqrt{1-x}}{x^2 \left (a+b x+c x^2\right )} \, dx,x,x^2\right )\\ &=-\operatorname{Subst}\left (\int \frac{x^2}{\left (1-x^2\right )^2 \left (a+b+c+(-b-2 c) x^2+c x^4\right )} \, dx,x,\sqrt{1-x^2}\right )\\ &=-\operatorname{Subst}\left (\int \left (\frac{1}{4 a (-1+x)^2}+\frac{1}{4 a (1+x)^2}+\frac{a+2 b}{2 a^2 \left (-1+x^2\right )}+\frac{b (a+b+c)-(a+b) c x^2}{a^2 \left (a+b+c-(b+2 c) x^2+c x^4\right )}\right ) \, dx,x,\sqrt{1-x^2}\right )\\ &=-\frac{1}{4 a \left (1-\sqrt{1-x^2}\right )}+\frac{1}{4 a \left (1+\sqrt{1-x^2}\right )}-\frac{\operatorname{Subst}\left (\int \frac{b (a+b+c)-(a+b) c x^2}{a+b+c+(-b-2 c) x^2+c x^4} \, dx,x,\sqrt{1-x^2}\right )}{a^2}-\frac{(a+2 b) \operatorname{Subst}\left (\int \frac{1}{-1+x^2} \, dx,x,\sqrt{1-x^2}\right )}{2 a^2}\\ &=-\frac{1}{4 a \left (1-\sqrt{1-x^2}\right )}+\frac{1}{4 a \left (1+\sqrt{1-x^2}\right )}+\frac{(a+2 b) \tanh ^{-1}\left (\sqrt{1-x^2}\right )}{2 a^2}+\frac{\left (c \left (a+b-\frac{b^2+a (b-2 c)}{\sqrt{b^2-4 a c}}\right )\right ) \operatorname{Subst}\left (\int \frac{1}{\frac{1}{2} (-b-2 c)-\frac{1}{2} \sqrt{b^2-4 a c}+c x^2} \, dx,x,\sqrt{1-x^2}\right )}{2 a^2}+\frac{\left (c \left (a+b+\frac{b^2+a (b-2 c)}{\sqrt{b^2-4 a c}}\right )\right ) \operatorname{Subst}\left (\int \frac{1}{\frac{1}{2} (-b-2 c)+\frac{1}{2} \sqrt{b^2-4 a c}+c x^2} \, dx,x,\sqrt{1-x^2}\right )}{2 a^2}\\ &=-\frac{1}{4 a \left (1-\sqrt{1-x^2}\right )}+\frac{1}{4 a \left (1+\sqrt{1-x^2}\right )}+\frac{(a+2 b) \tanh ^{-1}\left (\sqrt{1-x^2}\right )}{2 a^2}-\frac{\sqrt{c} \left (a+b+\frac{b^2+a (b-2 c)}{\sqrt{b^2-4 a c}}\right ) \tanh ^{-1}\left (\frac{\sqrt{2} \sqrt{c} \sqrt{1-x^2}}{\sqrt{b+2 c-\sqrt{b^2-4 a c}}}\right )}{\sqrt{2} a^2 \sqrt{b+2 c-\sqrt{b^2-4 a c}}}-\frac{\sqrt{c} \left (a+b-\frac{b^2+a (b-2 c)}{\sqrt{b^2-4 a c}}\right ) \tanh ^{-1}\left (\frac{\sqrt{2} \sqrt{c} \sqrt{1-x^2}}{\sqrt{b+2 c+\sqrt{b^2-4 a c}}}\right )}{\sqrt{2} a^2 \sqrt{b+2 c+\sqrt{b^2-4 a c}}}\\ \end{align*}

Mathematica [A]  time = 0.770169, size = 292, normalized size = 1.01 \[ \frac{\frac{\sqrt{2} \sqrt{c} \left (-\frac{\left (b \left (\sqrt{b^2-4 a c}+b\right )+a \left (\sqrt{b^2-4 a c}+b-2 c\right )\right ) \tanh ^{-1}\left (\frac{\sqrt{2} \sqrt{c} \sqrt{1-x^2}}{\sqrt{-\sqrt{b^2-4 a c}+b+2 c}}\right )}{\sqrt{-\sqrt{b^2-4 a c}+b+2 c}}-\frac{\left (b \left (\sqrt{b^2-4 a c}-b\right )+a \left (\sqrt{b^2-4 a c}-b+2 c\right )\right ) \tanh ^{-1}\left (\frac{\sqrt{2} \sqrt{c} \sqrt{1-x^2}}{\sqrt{\sqrt{b^2-4 a c}+b+2 c}}\right )}{\sqrt{\sqrt{b^2-4 a c}+b+2 c}}\right )}{\sqrt{b^2-4 a c}}+(a+2 b) \log \left (\sqrt{1-x^2}+1\right )-(a+2 b) \log (x)-\frac{a \sqrt{1-x^2}}{x^2}}{2 a^2} \]

Antiderivative was successfully verified.

[In]

Integrate[Sqrt[1 - x^2]/(x^3*(a + b*x^2 + c*x^4)),x]

[Out]

(-((a*Sqrt[1 - x^2])/x^2) + (Sqrt[2]*Sqrt[c]*(-(((b*(b + Sqrt[b^2 - 4*a*c]) + a*(b - 2*c + Sqrt[b^2 - 4*a*c]))
*ArcTanh[(Sqrt[2]*Sqrt[c]*Sqrt[1 - x^2])/Sqrt[b + 2*c - Sqrt[b^2 - 4*a*c]]])/Sqrt[b + 2*c - Sqrt[b^2 - 4*a*c]]
) - ((b*(-b + Sqrt[b^2 - 4*a*c]) + a*(-b + 2*c + Sqrt[b^2 - 4*a*c]))*ArcTanh[(Sqrt[2]*Sqrt[c]*Sqrt[1 - x^2])/S
qrt[b + 2*c + Sqrt[b^2 - 4*a*c]]])/Sqrt[b + 2*c + Sqrt[b^2 - 4*a*c]]))/Sqrt[b^2 - 4*a*c] - (a + 2*b)*Log[x] +
(a + 2*b)*Log[1 + Sqrt[1 - x^2]])/(2*a^2)

________________________________________________________________________________________

Maple [B]  time = 0.059, size = 2770, normalized size = 9.6 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int((-x^2+1)^(1/2)/x^3/(c*x^4+b*x^2+a),x)

[Out]

2/(4*a*c-b^2)/(4*a*c-2*b^2+2*(-4*a*c+b^2)^(1/2)*a+2*b*(-4*a*c+b^2)^(1/2)-2*a*b)^(1/2)*arctan(1/2*(-2*((-x^2+1)
^(1/2)-1)^2/x^2*a+2*(-4*a*c+b^2)^(1/2)-2*a-2*b)/(4*a*c-2*b^2+2*(-4*a*c+b^2)^(1/2)*a+2*b*(-4*a*c+b^2)^(1/2)-2*a
*b)^(1/2))*c*(-4*a*c+b^2)^(1/2)-1/a/(4*a*c-b^2)/(4*a*c-2*b^2+2*(-4*a*c+b^2)^(1/2)*a+2*b*(-4*a*c+b^2)^(1/2)-2*a
*b)^(1/2)*arctan(1/2*(-2*((-x^2+1)^(1/2)-1)^2/x^2*a+2*(-4*a*c+b^2)^(1/2)-2*a-2*b)/(4*a*c-2*b^2+2*(-4*a*c+b^2)^
(1/2)*a+2*b*(-4*a*c+b^2)^(1/2)-2*a*b)^(1/2))*b^2*(-4*a*c+b^2)^(1/2)+3/a/(4*a*c-b^2)/(4*a*c-2*b^2+2*(-4*a*c+b^2
)^(1/2)*a+2*b*(-4*a*c+b^2)^(1/2)-2*a*b)^(1/2)*arctan(1/2*(-2*((-x^2+1)^(1/2)-1)^2/x^2*a+2*(-4*a*c+b^2)^(1/2)-2
*a-2*b)/(4*a*c-2*b^2+2*(-4*a*c+b^2)^(1/2)*a+2*b*(-4*a*c+b^2)^(1/2)-2*a*b)^(1/2))*b*c*(-4*a*c+b^2)^(1/2)-1/a^2/
(4*a*c-b^2)/(4*a*c-2*b^2+2*(-4*a*c+b^2)^(1/2)*a+2*b*(-4*a*c+b^2)^(1/2)-2*a*b)^(1/2)*arctan(1/2*(-2*((-x^2+1)^(
1/2)-1)^2/x^2*a+2*(-4*a*c+b^2)^(1/2)-2*a-2*b)/(4*a*c-2*b^2+2*(-4*a*c+b^2)^(1/2)*a+2*b*(-4*a*c+b^2)^(1/2)-2*a*b
)^(1/2))*b^3*(-4*a*c+b^2)^(1/2)-4/(4*a*c-b^2)/(4*a*c-2*b^2+2*(-4*a*c+b^2)^(1/2)*a+2*b*(-4*a*c+b^2)^(1/2)-2*a*b
)^(1/2)*arctan(1/2*(-2*((-x^2+1)^(1/2)-1)^2/x^2*a+2*(-4*a*c+b^2)^(1/2)-2*a-2*b)/(4*a*c-2*b^2+2*(-4*a*c+b^2)^(1
/2)*a+2*b*(-4*a*c+b^2)^(1/2)-2*a*b)^(1/2))*b*c+4/(4*a*c-b^2)/(4*a*c-2*b^2+2*(-4*a*c+b^2)^(1/2)*a+2*b*(-4*a*c+b
^2)^(1/2)-2*a*b)^(1/2)*arctan(1/2*(-2*((-x^2+1)^(1/2)-1)^2/x^2*a+2*(-4*a*c+b^2)^(1/2)-2*a-2*b)/(4*a*c-2*b^2+2*
(-4*a*c+b^2)^(1/2)*a+2*b*(-4*a*c+b^2)^(1/2)-2*a*b)^(1/2))*c^2+1/a/(4*a*c-b^2)/(4*a*c-2*b^2+2*(-4*a*c+b^2)^(1/2
)*a+2*b*(-4*a*c+b^2)^(1/2)-2*a*b)^(1/2)*arctan(1/2*(-2*((-x^2+1)^(1/2)-1)^2/x^2*a+2*(-4*a*c+b^2)^(1/2)-2*a-2*b
)/(4*a*c-2*b^2+2*(-4*a*c+b^2)^(1/2)*a+2*b*(-4*a*c+b^2)^(1/2)-2*a*b)^(1/2))*b^3-5/a/(4*a*c-b^2)/(4*a*c-2*b^2+2*
(-4*a*c+b^2)^(1/2)*a+2*b*(-4*a*c+b^2)^(1/2)-2*a*b)^(1/2)*arctan(1/2*(-2*((-x^2+1)^(1/2)-1)^2/x^2*a+2*(-4*a*c+b
^2)^(1/2)-2*a-2*b)/(4*a*c-2*b^2+2*(-4*a*c+b^2)^(1/2)*a+2*b*(-4*a*c+b^2)^(1/2)-2*a*b)^(1/2))*b^2*c+1/a^2/(4*a*c
-b^2)/(4*a*c-2*b^2+2*(-4*a*c+b^2)^(1/2)*a+2*b*(-4*a*c+b^2)^(1/2)-2*a*b)^(1/2)*arctan(1/2*(-2*((-x^2+1)^(1/2)-1
)^2/x^2*a+2*(-4*a*c+b^2)^(1/2)-2*a-2*b)/(4*a*c-2*b^2+2*(-4*a*c+b^2)^(1/2)*a+2*b*(-4*a*c+b^2)^(1/2)-2*a*b)^(1/2
))*b^4+2/(4*a*c-b^2)/(4*a*c-2*b^2-2*(-4*a*c+b^2)^(1/2)*a-2*b*(-4*a*c+b^2)^(1/2)-2*a*b)^(1/2)*arctan(1/2*(2*((-
x^2+1)^(1/2)-1)^2/x^2*a+2*(-4*a*c+b^2)^(1/2)+2*a+2*b)/(4*a*c-2*b^2-2*(-4*a*c+b^2)^(1/2)*a-2*b*(-4*a*c+b^2)^(1/
2)-2*a*b)^(1/2))*c*(-4*a*c+b^2)^(1/2)-1/a/(4*a*c-b^2)/(4*a*c-2*b^2-2*(-4*a*c+b^2)^(1/2)*a-2*b*(-4*a*c+b^2)^(1/
2)-2*a*b)^(1/2)*arctan(1/2*(2*((-x^2+1)^(1/2)-1)^2/x^2*a+2*(-4*a*c+b^2)^(1/2)+2*a+2*b)/(4*a*c-2*b^2-2*(-4*a*c+
b^2)^(1/2)*a-2*b*(-4*a*c+b^2)^(1/2)-2*a*b)^(1/2))*b^2*(-4*a*c+b^2)^(1/2)+3/a/(4*a*c-b^2)/(4*a*c-2*b^2-2*(-4*a*
c+b^2)^(1/2)*a-2*b*(-4*a*c+b^2)^(1/2)-2*a*b)^(1/2)*arctan(1/2*(2*((-x^2+1)^(1/2)-1)^2/x^2*a+2*(-4*a*c+b^2)^(1/
2)+2*a+2*b)/(4*a*c-2*b^2-2*(-4*a*c+b^2)^(1/2)*a-2*b*(-4*a*c+b^2)^(1/2)-2*a*b)^(1/2))*b*c*(-4*a*c+b^2)^(1/2)-1/
a^2/(4*a*c-b^2)/(4*a*c-2*b^2-2*(-4*a*c+b^2)^(1/2)*a-2*b*(-4*a*c+b^2)^(1/2)-2*a*b)^(1/2)*arctan(1/2*(2*((-x^2+1
)^(1/2)-1)^2/x^2*a+2*(-4*a*c+b^2)^(1/2)+2*a+2*b)/(4*a*c-2*b^2-2*(-4*a*c+b^2)^(1/2)*a-2*b*(-4*a*c+b^2)^(1/2)-2*
a*b)^(1/2))*b^3*(-4*a*c+b^2)^(1/2)+4/(4*a*c-b^2)/(4*a*c-2*b^2-2*(-4*a*c+b^2)^(1/2)*a-2*b*(-4*a*c+b^2)^(1/2)-2*
a*b)^(1/2)*arctan(1/2*(2*((-x^2+1)^(1/2)-1)^2/x^2*a+2*(-4*a*c+b^2)^(1/2)+2*a+2*b)/(4*a*c-2*b^2-2*(-4*a*c+b^2)^
(1/2)*a-2*b*(-4*a*c+b^2)^(1/2)-2*a*b)^(1/2))*b*c-4/(4*a*c-b^2)/(4*a*c-2*b^2-2*(-4*a*c+b^2)^(1/2)*a-2*b*(-4*a*c
+b^2)^(1/2)-2*a*b)^(1/2)*arctan(1/2*(2*((-x^2+1)^(1/2)-1)^2/x^2*a+2*(-4*a*c+b^2)^(1/2)+2*a+2*b)/(4*a*c-2*b^2-2
*(-4*a*c+b^2)^(1/2)*a-2*b*(-4*a*c+b^2)^(1/2)-2*a*b)^(1/2))*c^2-1/a/(4*a*c-b^2)/(4*a*c-2*b^2-2*(-4*a*c+b^2)^(1/
2)*a-2*b*(-4*a*c+b^2)^(1/2)-2*a*b)^(1/2)*arctan(1/2*(2*((-x^2+1)^(1/2)-1)^2/x^2*a+2*(-4*a*c+b^2)^(1/2)+2*a+2*b
)/(4*a*c-2*b^2-2*(-4*a*c+b^2)^(1/2)*a-2*b*(-4*a*c+b^2)^(1/2)-2*a*b)^(1/2))*b^3+5/a/(4*a*c-b^2)/(4*a*c-2*b^2-2*
(-4*a*c+b^2)^(1/2)*a-2*b*(-4*a*c+b^2)^(1/2)-2*a*b)^(1/2)*arctan(1/2*(2*((-x^2+1)^(1/2)-1)^2/x^2*a+2*(-4*a*c+b^
2)^(1/2)+2*a+2*b)/(4*a*c-2*b^2-2*(-4*a*c+b^2)^(1/2)*a-2*b*(-4*a*c+b^2)^(1/2)-2*a*b)^(1/2))*b^2*c-1/a^2/(4*a*c-
b^2)/(4*a*c-2*b^2-2*(-4*a*c+b^2)^(1/2)*a-2*b*(-4*a*c+b^2)^(1/2)-2*a*b)^(1/2)*arctan(1/2*(2*((-x^2+1)^(1/2)-1)^
2/x^2*a+2*(-4*a*c+b^2)^(1/2)+2*a+2*b)/(4*a*c-2*b^2-2*(-4*a*c+b^2)^(1/2)*a-2*b*(-4*a*c+b^2)^(1/2)-2*a*b)^(1/2))
*b^4+2/a^2*b/(2/x^2-2/x^2*(-x^2+1)^(1/2))-1/a^2*b*(-x^2+1)^(1/2)+1/a^2*b*arctanh(1/(-x^2+1)^(1/2))-1/2/a/x^2*(
-x^2+1)^(3/2)-1/2/a*(-x^2+1)^(1/2)+1/2/a*arctanh(1/(-x^2+1)^(1/2))

________________________________________________________________________________________

Maxima [F]  time = 0., size = 0, normalized size = 0. \begin{align*} \int \frac{\sqrt{-x^{2} + 1}}{{\left (c x^{4} + b x^{2} + a\right )} x^{3}}\,{d x} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((-x^2+1)^(1/2)/x^3/(c*x^4+b*x^2+a),x, algorithm="maxima")

[Out]

integrate(sqrt(-x^2 + 1)/((c*x^4 + b*x^2 + a)*x^3), x)

________________________________________________________________________________________

Fricas [B]  time = 67.3797, size = 5650, normalized size = 19.48 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((-x^2+1)^(1/2)/x^3/(c*x^4+b*x^2+a),x, algorithm="fricas")

[Out]

-1/2*(sqrt(1/2)*a^2*x^2*sqrt((a*b^3 + b^4 + 2*a^2*c^2 - (3*a^2*b + 4*a*b^2)*c - (a^4*b^2 - 4*a^5*c)*sqrt((a^2*
b^4 + 2*a*b^5 + b^6 + (a^4 + 4*a^3*b + 4*a^2*b^2)*c^2 - 2*(a^3*b^2 + 3*a^2*b^3 + 2*a*b^4)*c)/(a^8*b^2 - 4*a^9*
c)))/(a^4*b^2 - 4*a^5*c))*log(((a^4*b^2*c - 4*a^5*c^2)*x^2*sqrt((a^2*b^4 + 2*a*b^5 + b^6 + (a^4 + 4*a^3*b + 4*
a^2*b^2)*c^2 - 2*(a^3*b^2 + 3*a^2*b^3 + 2*a*b^4)*c)/(a^8*b^2 - 4*a^9*c)) + 2*(a^3 + 2*a^2*b)*c^2 + ((a^2*b + 2
*a*b^2)*c^2 - (a*b^3 + b^4)*c)*x^2 - 2*(a^2*b^2 + a*b^3)*c + sqrt(1/2)*((a^5*b^3 - 4*a^6*b*c)*x^2*sqrt((a^2*b^
4 + 2*a*b^5 + b^6 + (a^4 + 4*a^3*b + 4*a^2*b^2)*c^2 - 2*(a^3*b^2 + 3*a^2*b^3 + 2*a*b^4)*c)/(a^8*b^2 - 4*a^9*c)
) + (a^2*b^4 + a*b^5 + 4*(a^4 + 2*a^3*b)*c^2 - (5*a^3*b^2 + 6*a^2*b^3)*c)*x^2)*sqrt((a*b^3 + b^4 + 2*a^2*c^2 -
 (3*a^2*b + 4*a*b^2)*c - (a^4*b^2 - 4*a^5*c)*sqrt((a^2*b^4 + 2*a*b^5 + b^6 + (a^4 + 4*a^3*b + 4*a^2*b^2)*c^2 -
 2*(a^3*b^2 + 3*a^2*b^3 + 2*a*b^4)*c)/(a^8*b^2 - 4*a^9*c)))/(a^4*b^2 - 4*a^5*c)) - 2*((a^3 + 2*a^2*b)*c^2 - (a
^2*b^2 + a*b^3)*c)*sqrt(-x^2 + 1))/x^2) - sqrt(1/2)*a^2*x^2*sqrt((a*b^3 + b^4 + 2*a^2*c^2 - (3*a^2*b + 4*a*b^2
)*c - (a^4*b^2 - 4*a^5*c)*sqrt((a^2*b^4 + 2*a*b^5 + b^6 + (a^4 + 4*a^3*b + 4*a^2*b^2)*c^2 - 2*(a^3*b^2 + 3*a^2
*b^3 + 2*a*b^4)*c)/(a^8*b^2 - 4*a^9*c)))/(a^4*b^2 - 4*a^5*c))*log(((a^4*b^2*c - 4*a^5*c^2)*x^2*sqrt((a^2*b^4 +
 2*a*b^5 + b^6 + (a^4 + 4*a^3*b + 4*a^2*b^2)*c^2 - 2*(a^3*b^2 + 3*a^2*b^3 + 2*a*b^4)*c)/(a^8*b^2 - 4*a^9*c)) +
 2*(a^3 + 2*a^2*b)*c^2 + ((a^2*b + 2*a*b^2)*c^2 - (a*b^3 + b^4)*c)*x^2 - 2*(a^2*b^2 + a*b^3)*c - sqrt(1/2)*((a
^5*b^3 - 4*a^6*b*c)*x^2*sqrt((a^2*b^4 + 2*a*b^5 + b^6 + (a^4 + 4*a^3*b + 4*a^2*b^2)*c^2 - 2*(a^3*b^2 + 3*a^2*b
^3 + 2*a*b^4)*c)/(a^8*b^2 - 4*a^9*c)) + (a^2*b^4 + a*b^5 + 4*(a^4 + 2*a^3*b)*c^2 - (5*a^3*b^2 + 6*a^2*b^3)*c)*
x^2)*sqrt((a*b^3 + b^4 + 2*a^2*c^2 - (3*a^2*b + 4*a*b^2)*c - (a^4*b^2 - 4*a^5*c)*sqrt((a^2*b^4 + 2*a*b^5 + b^6
 + (a^4 + 4*a^3*b + 4*a^2*b^2)*c^2 - 2*(a^3*b^2 + 3*a^2*b^3 + 2*a*b^4)*c)/(a^8*b^2 - 4*a^9*c)))/(a^4*b^2 - 4*a
^5*c)) - 2*((a^3 + 2*a^2*b)*c^2 - (a^2*b^2 + a*b^3)*c)*sqrt(-x^2 + 1))/x^2) + sqrt(1/2)*a^2*x^2*sqrt((a*b^3 +
b^4 + 2*a^2*c^2 - (3*a^2*b + 4*a*b^2)*c + (a^4*b^2 - 4*a^5*c)*sqrt((a^2*b^4 + 2*a*b^5 + b^6 + (a^4 + 4*a^3*b +
 4*a^2*b^2)*c^2 - 2*(a^3*b^2 + 3*a^2*b^3 + 2*a*b^4)*c)/(a^8*b^2 - 4*a^9*c)))/(a^4*b^2 - 4*a^5*c))*log(-((a^4*b
^2*c - 4*a^5*c^2)*x^2*sqrt((a^2*b^4 + 2*a*b^5 + b^6 + (a^4 + 4*a^3*b + 4*a^2*b^2)*c^2 - 2*(a^3*b^2 + 3*a^2*b^3
 + 2*a*b^4)*c)/(a^8*b^2 - 4*a^9*c)) - 2*(a^3 + 2*a^2*b)*c^2 - ((a^2*b + 2*a*b^2)*c^2 - (a*b^3 + b^4)*c)*x^2 +
2*(a^2*b^2 + a*b^3)*c + sqrt(1/2)*((a^5*b^3 - 4*a^6*b*c)*x^2*sqrt((a^2*b^4 + 2*a*b^5 + b^6 + (a^4 + 4*a^3*b +
4*a^2*b^2)*c^2 - 2*(a^3*b^2 + 3*a^2*b^3 + 2*a*b^4)*c)/(a^8*b^2 - 4*a^9*c)) - (a^2*b^4 + a*b^5 + 4*(a^4 + 2*a^3
*b)*c^2 - (5*a^3*b^2 + 6*a^2*b^3)*c)*x^2)*sqrt((a*b^3 + b^4 + 2*a^2*c^2 - (3*a^2*b + 4*a*b^2)*c + (a^4*b^2 - 4
*a^5*c)*sqrt((a^2*b^4 + 2*a*b^5 + b^6 + (a^4 + 4*a^3*b + 4*a^2*b^2)*c^2 - 2*(a^3*b^2 + 3*a^2*b^3 + 2*a*b^4)*c)
/(a^8*b^2 - 4*a^9*c)))/(a^4*b^2 - 4*a^5*c)) + 2*((a^3 + 2*a^2*b)*c^2 - (a^2*b^2 + a*b^3)*c)*sqrt(-x^2 + 1))/x^
2) - sqrt(1/2)*a^2*x^2*sqrt((a*b^3 + b^4 + 2*a^2*c^2 - (3*a^2*b + 4*a*b^2)*c + (a^4*b^2 - 4*a^5*c)*sqrt((a^2*b
^4 + 2*a*b^5 + b^6 + (a^4 + 4*a^3*b + 4*a^2*b^2)*c^2 - 2*(a^3*b^2 + 3*a^2*b^3 + 2*a*b^4)*c)/(a^8*b^2 - 4*a^9*c
)))/(a^4*b^2 - 4*a^5*c))*log(-((a^4*b^2*c - 4*a^5*c^2)*x^2*sqrt((a^2*b^4 + 2*a*b^5 + b^6 + (a^4 + 4*a^3*b + 4*
a^2*b^2)*c^2 - 2*(a^3*b^2 + 3*a^2*b^3 + 2*a*b^4)*c)/(a^8*b^2 - 4*a^9*c)) - 2*(a^3 + 2*a^2*b)*c^2 - ((a^2*b + 2
*a*b^2)*c^2 - (a*b^3 + b^4)*c)*x^2 + 2*(a^2*b^2 + a*b^3)*c - sqrt(1/2)*((a^5*b^3 - 4*a^6*b*c)*x^2*sqrt((a^2*b^
4 + 2*a*b^5 + b^6 + (a^4 + 4*a^3*b + 4*a^2*b^2)*c^2 - 2*(a^3*b^2 + 3*a^2*b^3 + 2*a*b^4)*c)/(a^8*b^2 - 4*a^9*c)
) - (a^2*b^4 + a*b^5 + 4*(a^4 + 2*a^3*b)*c^2 - (5*a^3*b^2 + 6*a^2*b^3)*c)*x^2)*sqrt((a*b^3 + b^4 + 2*a^2*c^2 -
 (3*a^2*b + 4*a*b^2)*c + (a^4*b^2 - 4*a^5*c)*sqrt((a^2*b^4 + 2*a*b^5 + b^6 + (a^4 + 4*a^3*b + 4*a^2*b^2)*c^2 -
 2*(a^3*b^2 + 3*a^2*b^3 + 2*a*b^4)*c)/(a^8*b^2 - 4*a^9*c)))/(a^4*b^2 - 4*a^5*c)) + 2*((a^3 + 2*a^2*b)*c^2 - (a
^2*b^2 + a*b^3)*c)*sqrt(-x^2 + 1))/x^2) + (a + 2*b)*x^2*log((sqrt(-x^2 + 1) - 1)/x) + sqrt(-x^2 + 1)*a)/(a^2*x
^2)

________________________________________________________________________________________

Sympy [F(-1)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Timed out} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((-x**2+1)**(1/2)/x**3/(c*x**4+b*x**2+a),x)

[Out]

Timed out

________________________________________________________________________________________

Giac [F(-1)]  time = 0., size = 0, normalized size = 0. \begin{align*} \text{Timed out} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate((-x^2+1)^(1/2)/x^3/(c*x^4+b*x^2+a),x, algorithm="giac")

[Out]

Timed out